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FOREWORD 


This  note  describes  the  results  of  an  attempt  to  compute  the  envelope  and 
frequency  of  an  analytic  signal  from  a set  of  sampled  data.  This  note  has  been 
prepared  since  it  may  be  of  interest  to  a limited  number  of  persons  at  the  U.  S. 
Navy  Undersea  Warfare  Center.  It  should  not  be  constriied  as  a formal  report  since 
its  function  is  to  present  for  information  a small  portion  of  work  which  is  ap- 
plicable to  the  area  of  sonar  signal  measurements  and  analysis. 

This  note  was  written  as  a term  paper  for  a course  in  signal  analysis  at 
Purdue  University  in  the  Fall  of  19 66.  The  paper  presented  the  results  of  an 
analysis  of  an  interferogram  waveform.  It  was  assumed  at  that  time  that  the 
reader  had  intimate  acquaintance  with  such  a device.  Therefore,  the  method  of 
signal  generation  was  not  discussed.  In  making  this  paper  available  as  a NUWC 
Technical  note,  it  is  now  necessary  to  include  a short  description  of  the  inter- 
ferogram signal  in  Appendix  A. 
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1.  Introduction. 


This  report  presents  the  results  of  an  attempt  to  compute  the 
envelope  and  frequency  of  an  analytic  signal  from  a set  of  sampled 
data.  Interferogram  data  from  run  4 (flashlight  bulb)  were  used  as 
a data  base.  Computations  were  performed  on  the  Purdue  IBM  7094. 

2.  Background. 

2.1  The  analytic  signal. 

The  definition  of  envelope  and  frequency  may  be  obtained  by 
considering  the  analytic  signal  of  the  form 

/*■ 

f(t)  = v (t)  e (2-i) 

The  envelope  is  then  the  real  time  varying  function  V(t)  where  V(t)>0. 

is  a complex  function,  a phasor,  which  rotates  with  angular 
(radian)  frequency  4>(t). 

The  question  is,  if  we  are  presented  with  an  analog  signal,  f(t), 
which  is  a real  function  of  time,  how  can  we  obtain  envelope  and  frequency 
information  from  it? 

The  procedure  used  here  is  to  consider  another  form  of  the  analytic 
signal  which  is  equivalent  to  the  one  given  above.  That  is, 

f(t)  = f(t)  + j f(t) 


(2-2) 


where  f(t)  is  the  Hilbert  transform  of  f(t)  (written  as  H{f(t)}  and 
called  the  quadrature  function  of  f(t))  and  ll{f(t)}  is  defined  as* 

«<'(')>  (2-5’ 

'-CO 

The  equivalence  is  shown  by  the  following  relations: 


V(t)  = (f*(t)  ♦ f‘(t))''1 


and 


4>(t) 

*(t) 


tan’ 

Jt 


A 

ft*) 


2.2  The  sampling  theorem 


(2-4) 

(2-5) 

(2-6) 


The  sampling  theorem  states  that,  given  a time  function,  f(t), 
which  contains  no  frequencies  higher  than  IV  Hz,  and  also  given  a set 
of  samples  of  the  function  spaced  no  more  than  1/2  W seconds  apart,  it 
is  possible  to  recover  the  intervening  values  of  the  time  function  with 
full  accuracy.  Reconstruction  of  the  original  time  function,  f(t),  can 
be  accomplished  by  using  the  following  relation: 


f(t)  «=  2_ f(nAt) 


jiu  (tt-^  - /w^f) 
(it  £ -Aor) 


(2-7) 


*Onc  has  to  be  careful  when  referring  to  some  texts  for  tables  of 
Hilbert  transforms.  There  are  a number  of  ways  of  defining  the 
Hilbert  transform,  all  depending  on  whether  one  wants  to  end  up 
with  the  original  function  or  its  negative  after  two  Hilbert 
transformat  ions  in  succession. 
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where  At  *=  1/2  W is  the  time  interval  between  samples  and  n «=  ...-1,  0, 


♦1. 


To  see  why  this  relation  works,  we  will  first  look  at  the  spectrum 
of  the  sampled  signal.  In  effect,  the  sampling  process  may  be  considered 
as  equivalent  to  multiplying  f(t)  by  a train  of  delta- functions  spaced 
At  seconds  apart.  The  Fourier  transform  of  such  a product  is  a train 
of  delta- functions  in  the  frequency  domain  convolved  with  the  spectrum 
of  the  original  time  function  f(t)*.  That  is, 

go 

5>  - nAt)f(t)  ♦ m ~ ) * F(w) 

where 

f(nAt)  = f(t)  g(t  - nAt)  (2-S) 

and 

f(t)  <i=£>  F(w). 

The  frequency  spectrum  of  f(nAt)  consists  of  a train  of  replicas  of 
F(u>)  spaced  at  positive  and  negative  integer  multiples  of  QK/tt  radians. 
Therefore,  in  order  to  recover  the  original  spectrum  (and  assuming  spec- 
trum overlap  between  adjacent  spectrum  replicas  doesn't  occur)**  it  is 

*A.  Papoulis,  The  Fourier  Integral  and  its  Application,  McGraw-Hill, 
1962,  p.  48. 

**This  is  the  equivalent  to  the  restriction  given  in  the  statement  of  the 
sampling  theorem  where  the  samples  must  be  spaced  no  more  than  1/2  w 
seconds  apart. 
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sufficient  to  pass  the  spectrum  centered  at  u>  = 0 and  reject  all  other 
replicas  by  using  an  ideal  rectangular  band  pass  filter  with  the 
characteristics 


H(w) 


1 |wl<6 

O elsewhere. 


The  impulse  response  of  this  filter  is  simply 

Sin  0tr 


1(0 


0t 


(2-9) 


Thus,  when  the  time  function  f(nAt)  is  passed  through  the  filter 
the  output  waveform  is  the  convolution  of  f(nAt)  with  the  impulse 
response  of  the  filter.  Then 


f(t) 


SiM  Bt' 
iit 


= J”f(nAt)  ’ 

V'  , sikj  6(t  - 
= 1 f(nAt)  — 


(2-10) 

(2-11) 


ft  (t  - >h/3W) 

If  we  set  B = 2HW  = then  equation  (2-11)  becomes  equation  (2-7). 
The  relation  B *=  2IW  is  a convenient  one  since  this  point  is  midway 
between  the  center  of  P(w)  and  its  adjacent  replica. 

3.  Computation  of  envelope  and  frequency. 

Equation  (2-7)  gives  us  a relation  between  f(t)  (the  reconstructed 
function)  and  the  sampled  values  of  f(t)  (the  original  function),  f(nAt) 

A 

What  we  now  need  is  the  quadrature  function  f(t).  This  can  be  obtained 
as  follows: 
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f(t)  = H{f (t) } 


(3-1) 


( 5im  ( if  & - /HTf  ) 

H ) f(nAt) 


n 


(IT  - /W7T 

Hf 

fe-  I (n't  - /mtt  ) J 

At  ' 

I - CoS  (if  SI  "A**  ) 


(3-2) 


(3-3) 


= ^~f (nAt) 


(tt^  -a.it) 

The  sin(x)/x  function  and  its  Hilbert  transform  are  shown  in 
Figure  1 for  comparison. 


(3-4) 


The  envelope  and  phase  are  computed  using  equations  (2-4),  (2-5), 
(2-7),  and  (3-4). 


Rad  ian  frequency  is  computed  by  using  the  approximation 


♦ (t  • ) * 


4>te)  - ) 


*-•  - tj., 

where  t;  and  t^.,  arc  two  neighboring  points  in  time  for  which  f(t) 


and  f(t)  are  computed. 


4 . Results  of  computer  analysis . 


Figures  2 and  3 show  the  results  of  the  reconstruction  over  the 
interval  (sample  number)  393  seconds  to  435  seconds.  This  interval 
was  sufficient  to  enclose  the  central  peak  of  the  ’nterferogram  wave- 
form. All  waveforms  in  the  figures  were  computed  using  original  data 
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samples  numbered  375  through  455.  In  other  words,  in  equations  (2-11) 
and  (3-4)  n was  varied  from  375  to  455  while  t was  varied  from  393  to 
435.  Rough  indications  are  that  the  accuracy  of  the  computed  amplitudes 
is  within  2%  and  the  accuracy  of  the  computed  frequency  is  within  1%  of 
the  values  that  would  result  if  n varied  over  a much  wider  range  (say 
for  n<355  and  nM75)  . 

Figures  4 and  5 show  the  results  of  reconstruction  over  the  interval 
230  to  250  seconds.  In  this  case  n was  varied  from  200  to  640  and  t was 
varied  from  230  to  250. 


The  interferogram  waveform  appears  to  have  remained  negative  in 
value  over  most  of  the  interval.  This  is  due  to  the  12  Hz  noise  which 
was  added  to  the  interferometer  waveform.  Also  the  computed  frequency 
is  a peculiar  shape.  Note  that  the  positive  peaks  in  frequency  are 
rounded  and  correspond  to  the  greatest  (most  negative)  excursions  of 
the  time  function  or  the  greatest  (most  positive)  excursions  of  its 
Hilbert  transform,  whichever  happens  to  predominate.  Whereas,  the  nega- 
tive peaks  are  sharper  and  more  jagged  than  the  positive  peaks.  They 
also  correspond  in  time  to  the  lesser  (least  negative)  excursions  of 
the  time  function  or  the  lesser  (least  positive)  excursions  of  its 
Hilbert  transform,  whichever  happens  to  predominate. 

These  effects  can  be  explained  if  we  consider  the  following  wave- 
form: 
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' If 


f(t)  = A, cos  io,t  + AjCos  iD„t 


and  its  Hilbert  transform 

A 

f(t)  = A, sin  id,  t + A, sin  tD,t 

then 


f(t)  c A,{cos  u.t  + j sin  ui,t) 

♦ At{cos  w4t  + j sin  »tt) 

If  we  set  A,>>A4  and  u,<<w4  we  have  a situation  which  corresponds  quite 
closely  with  that  existing  with  the  interferogram  data.  Figure  6 
shows  a phasor  representation  for  the  above  equations,  h'e  sec  that 

A* 

maximum  velocity  of  a point  on  the  tip  of  the  phasor  representing  f(t) 
takes  place  when  A,  and  A4  are  colinear  and  pointing  in  the  same 
direction.  Conversely,  minimum  velocity  takes  place  when  A , and  A4  are 
colinear  but  this  time  pointing  in  opposite  directions.  This  minimum 
velocity  can  actually  be  negative  if  A,w,>  A4w4.  The  maximum  (minimum) 
velocity  takes  place  when  |f(t)(  is  a maximum  (minimum)  and  f(t)  is  zero 
when  A,  is  pointing  along  the  real  axis.  Conversely  the  maximum  (minimum) 
velocity  takes  place  when  |f(t)|  is  a maximum  (minimum)  and  f(t)  is  zero 
when  A,  is  pointing  along  the  imaginary  axis. 


From  Figure  6 we  can  see  that  the  maximum  and  minimum  for  the 
frequency  of  f(t)  is 
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maximum  frequency:  4>(t)  = 


4,  Ui  t AiliJi 


A i + Ai 


minimum  frequency: 


5(0 


A,  cO,  - 


A,  - Ax 


Since  A,,  Aj,  to,,  ut>  0 the  maximum  frequency  will  always  remain 
well  behaved.  However,  for  the  minimum  frequency,  when  A,  is  approximately 
equal  to  (A t - t ) the  denominator  becomes  arbitrarily  small.  ( e is  a 
small  number  greater  than  0).  If  u,<  we  observe  a large  negative- 
going "spike"  in  the  frequency. 

The  overall  appearance  of  the  waveform  f(t)  is  that  of  a coil 
rotating  counterclockwise  at  a rate  of  u, radians  per  second.  At  the 
tip  of  the  phasor  describing  this  coil  is  another  phasor  tracing  out 
another  smaller  coil  at  a faster  rate  of  radians  per  second. 

This  coil  too  is  rotating  in  a counterclockwise  direction.  This 
description  also  appears  to  be  appropriate  for  the  interferogram  data 
in  Figures  4 and  5. 

5.  Discussion. 

Until  now  the  subject  treated  here  was  primarily  an  outline  of  a 
technique  for  computing  the  envelope  and  frequency  of  an  analytic 
signal  as  based  on  a set  of  sampled  data,  and  an  exposure  to  some  of 
the  experimental  results  that  can  be  obtained  using  the  technique. 

Hardly  touched  on  was  the  subject  of  exceptions  to  some  of  the 
hypotheses  used  for  generating  the  technique.  Also  glossed  over  was 
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mention  of  errors  and  their  sources.  And  last  but  not  least  is  the 
question  of  the  utility  of  the  analytic  signal.  What  we  shall  do  at 
this  point  is  to  list  some  of  the  above  subjects  especially  as  they 
can  be  related  to  the  data  at  and.  What  we  will  end  up  with  is  a 
series  of  questions  with  partial  and  incomplete  answers  which  suggest 
areas  for  future  investigation. 

5. 1 Of  what  utility  is  the  analytic  signal  computation  technique? 

Perhaps  it  is  best  to  talk  about  analytic  waveforms  only  for  narrow- 
band  signals,  and  then  only  when  a specific  application  exis’s.  For 
example,  it  is  not  known  what  utility  exists  for  interferogra n wave- 
forms. Examples  can  be  drawai  from  other  areas,  however.  Sonar  and 
radar  returns  (echoes)  many  times  have  components  of  doppler  associated 
with  them.  Perhaps  this  method  of  computing  frequency  would  be  useful. 

5 . 2 What  are  the  exceptions  to  the  hypotheses  used  for  the  sampling 
theorem?  The  primary  exception  is  the  requirement  for  a band  limited 
signal.  Band  limited  signals  do  not  exist  in  nature.  In  general,  physical 
signals  exist  only  over  a finite  interval.  The  spectra  of  such  signals 
contain  frequency  components  at  all  frequencies.  However,  a spectrum 

may  contain  energy  concentrated  at  some  central  frequency  so  as  to  be 
approximately  band  limited.  In  such  a case  the  energy  outside  the 
band  limits  may  be  "negligible"  or  the  signal  may  be  filtered  to  cut 
this  energy  down  to  a negligible  level.  There  still  remains  the 
problem  of  estimating  the  amount  of  variation  in  the  computation  of 
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envelope  and  frequency  of  such  data  as  a function  of  the  amount  of 
spectrum  falling  outside  the  band  limits.  In  other  words,  how 
negligible  is  a good  question.  How  much  accuracy  do  we  sacrifice  using 
this  assumption? 

5.3  What  are  the  effects  of  data  acquisition  devices  on  the  accuracy 
of  the  computed  envelope  and  frequency?  This  is  impossible  to  estimate 
for  the  interferogram  data.  But  we  can  at  least  point  out  some  of  the 
sources  of  error.  Finite  quantization  levels,  finite  width  aperture 
of  the  sampler,  jitter  in  the  sampling  time,  and  circuit  noise  all 
contribute  to  errors  in  measurement.  Lumped  together  the  total  effect 
is  equal  to  that  of  an  equivalent  noise  source  in  series  with  an  ideal 
noiseless  signal.  We  also  assume  that  an  ideal  infinite  resolution 
sampler  is  used  for  measuring  the  signal  plus  noise.  Thus,  when 
reconstruction  is  attempted,  the  reconstructed  value  of  the  functions 

A 

f(t)  and  f(t)  will  each  have  errors  contributed  to  them  from  all  the 
time  series  samples  used  in  the  reconstruction.  These  errors  are 
carried  through  in  the  computation  of  envelope  and  frequency.  We  can 
say  that  those  time  series  samples  closest  to  the  reconstructed 
functions  in  time  will  contribute  not  only  the  greatest  amount  to 
reconstruction,  but  also  the  greatest  amount  to  the  error.  Indications 
are  that  the  reconstructed  value  of  the  functions  are  subject  to  the 
same  errors  as  each  of  the  sampled  data  points.*  This  subject  needs 

*R.  Braccwell,  The  Fourier  Transform  and  its  Applications,  McGraw-Hill, 
1965,  p.  203. 
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further  investigation.  Regardless  of  the  amount  of  added  noise  the 
computed  envelope  and  frequency  will  appear  to  have  a smooth  appearance. 
Figures  4 and  5 are  good  examples.  The  quantization  levels  of  the 
sampled  data  are  quite  evident  on  the  graph  of  the  computed  time 
function.  Yet  the  reconstruction  is  quite  smooth. 

5 . 4 What  effect  do  finite  data  samples  have  on  the  accuracy  of 
reconstruction?  The  problem  is  this.  Say  we  arc  interested  in 
reconstructing  the  envelope  and  frequency  of  a waveform  over  a time 
interval  T seconds  long.  Contributions  to  the  reconstruction  of 
the  function  at  any  point  in  time  within  that  interval  is  made  by 
neighboring  time  series  samples  on  both  sides  of  the  reconstruction 
point.  The  closer  in  time  the  samples  are  to  the  reconstruction 
point  the  more  they  contribute  to  the  value  of  the  function.*  Let 
us  now  bracket  the  reconst l'uct ion  point  by  a time  interval  T'  seconds 
long.  Then  suppose  we  perfoi'in  the  reconstruction  using  time  series 
samples  within  seconds  of  the  reconstruction  point.  How  large  must 

T*  be  in  order  to  reconstruct  f(t)  and  f(t)  to  within  a specified 
accuracy?  Two  reasons  for  answering  this  question  can  be  given.  The 
first  is  that  if  our  accuracy  is  going  to  be  degraded  some  certain 
amount  due  to  the  effects  of  the  equivalent  noise  discussed  in  part 

*Thi s is  true  in  general  for  both  f(t)  and  f(t).  However,  if  the 
reconstruction  point  is  the  same  point  in  time  as  a time  scries 
sample,  contributions  to  f(t)  come  only  from  that  time  series  sample 
and  no  other.  Whereas,  contributions  to  ?(t.)  come  from  all  time 
series  samples  except  that  one.  See  Figure  1. 
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5.  There  is  not  much  point  in  carrying  out  computations  to  a much 
greater  accuracy  by  including  too  much  data.  Also  computations  of 
the  sort  necessary  for  reconstruction  take  an  appreciable  amount  of 
computer  time,  and  time  costs  money.  Therefore  it  would  be  an 
advantage  to  be  able  to  estimate  T;  . Section  4 discussed  some 
results  where  the  most  accurate  reconstruction  of  the  functions  f(t) 

A 

and  f(t)  within  the  interval  395  to  435  seconds  was  obtained  from 
data  within  the  interval  355  to  475  seconds,  hooking  at  the  end 
points  of  the  smaller  interval  we  find  T'  is  equal  to  about  40  seconds 
(or  data  points  since  the  points  were  set  at  one  sample  per  second). 

5.5  How  much  error  is  introduced  by  the  digital  computer  computations? 
Wc  have  no  idea  how  much  error  is  introduced  but  we  can  point  to 
round-off  error  as  being  possibly  the  chief  source  of  error. 

6.  Conclusions 
\ ~ 

Computation  of  envelope  and  frequency  of  an  analytic  signal  based 
on  sampled  data  has  been  shown  to  be  feasible.  However,  there  exist 
a number  of  problem  areas  that  need  investigation  before  wc  can 
specify  the  amount  of  accuracy  that  these  computations  yield.  These 
areas  constitute  a list  of  possible  sources  of  error.* 

r*-  tfhc  error  due  to  the  use  of  non-band  limited  signals. 

- / 

b.  fhe  error  caused  by  the  equivalent  noise  associated  with  data 
acquisition  devices.  — ' • . « 
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FIGURE  1.  Comparison  of  - ' an<3  its  Hilbert  transform. 


Sample  number 


FIGURE  3.  Computed  Hilbert  transform  of  the  computed  interferogram  waveform  shown  in  Figure  2 


-V(t)  amplitude  of  computed  envelope 


I 


Instantaneous  frequency  r <j>(t) 

. radial  velocity  of  resultant 
radius  of  resultant 


FIGURE  6.  Phaser  representation  of  the  analytic  form  of  the 
sum  of  two  sine  waves. 
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APPENDIX  A 


A SHORT  DESCRIPTION  OF  THE  INTF.RFEROGRAM  SIGNAL 

A Michelson  interferometer  was  used  to  generate  an  interferogram 
signal  which  was  then  recorded  on  the  Ampex  SP  300  analog  tape  recorder. 
Also  recorded  was  a clock  waveform,  a 2 kHz  sine  wave,  which  was  used 
on  playback  to  trigger  an  analog-to-digital  converter  for  time  series 
samples  of  the  replayed  interferogram  waveform.  For  the  purpose  of 
the  analysis  it  was  assumed  that  the  sample  rate  was  one  per  second, 
not  the  actual  two  thousand  samples  per  second.  (Only  the  envelope 
and  frequency  computation  technique  was  of  interest  here.  The  inter- 
ferogram waveform  was  just  a convenient  source  of  data.)  About  950 
samples  were  taken  of  the  one  waveform  used  in  this  paper.  Not  all  of 
them  were  used  in  the  analysis.  Figure  A-l  illustrates  the  general 
characteristics  of  this  waveform. 


